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ABSTRACT 

ISO observed, for the first time to such a high sensitiv¬ 
ity level, the mid- and far- infrared universe. A Number 
of deep surveys were performed to probe the cosmological 
evolution of galaxies. In this review, I discuss and sum¬ 
marize results of mid-infrared ISOCAM and far-infrared 
ISOPHOT surveys, and show how our vision of the extra- 
galactic infrared universe has become more accurate. In 
particular, ISO allowed us to resolve into sources a signif¬ 
icant fraction of the Cosmic Infrared Background (CIB) 
in the mid-infrared, and to probe a fainter population in 
the far-infrared with the detection of the CIB fluctuations. 
Together with other wavelength data sets, the nature of 
ISO galaxies is now in the process of being understood. 

I also show that the high quality of the ISO data put 
strong constraints on the scenarios of galaxy evolution. 
This induced a burst in the development of models, yield¬ 
ing to a more coherent picture of galaxy evolution. 

I finally emphasize the potential of the ISO data archive 
in the field of observational cosmology, and describe the 
next steps, in particular the forthcoming cosmological sur¬ 
veys to be carried out by SIRTF. 
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extragalactic surveys - cosmic infrared background 


1. INTRODUCTION 

The Cosmic Infrared Background (CIB) (e.g. Puget et al, 
1996; Hauser & Dwek, 2001) is the relic emission of galaxy 
formation and evolution, being composed of light radiated 
from galaxies since their formation. The aim of the cos¬ 
mological surveys is to resolve the CIB into sources, and 
to provide data on galaxies at various wavelengths, red- 
shifts, with enough quality on individual objects as well 
as with statistical significant samples, in order to address 
the questions of the processes of galaxy formation and 
evolution, and the nature of the galaxies. 

The CIB Spectral Energy Distribution (SED), shown 
in Fig. 1 (from Hauser & Dwek, 2001) shows the existence 
of a minimum between 3 and 10 /rm separating direct stel¬ 
lar radiation from the FIR part due to radiation re-emitted 
by dust. The latter radiation contains at least a compara¬ 
ble integrated power as the former, and perhaps as much 
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Figure 1. The Cosmic Infrared Background Spectral En¬ 
ergy Distribution from Hauser & Dwek (2001). This plot 
includes detection (solid line in the submm part), upper 
limits (triangles and diamonds), and lower limits from 
source counts (sguares). The shaded region represents cur¬ 
rent observational limits, and the solid rising curve in the 
mm range represents the CMB.. 


as 2.5 times more. This ratio is much larger than what 
is measured locally (^30%). The CIB is thus likely to be 
dominated by a population of strongly evolving redshifted 
IR galaxies. 

ISO revolutionized our view of the extragalactic uni¬ 
verse with sensitive cosmological surveys carried out in 
the mid- and far- infrared. Of course, other cosmological 
surveys have been carried out with SCUBA, MAMBO, 
Chandra, XMM, and in the optical range; all give invalu¬ 
able informations about galaxy evolution, but it is beyond 
the scope of this short paper to review them. Reviews 
of extragalactic results from ISO can be found in Genzel 
& Cesarsky (2000) and Franceschini et al. (2001), and a 
comprehensive overview of ISO operations and science in 
Casoli et al. (2000) (particularly in the Puget & Guider- 
doni paper). I will here summarize the main results ob¬ 
tained in the MIR with ISOCAM (Sect. 2), in the FIR 
with ISOPHOT (Sect. 3), and overview the models that 
ISO data were able to constrain (Sect. 4). Sect. 5 empha¬ 
sizes the potential of the ISO archive, and in Sect. 6 I 
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Figure 2. K-corrections from Francecshini et al (2001). 
Dash: M82 spectrum at 7 /xm. Dot: inactive spiral at 15 
p,m. Line: M82-like spectrum at 15 fj,m. 

briefly explore the future of IR (and mm) astronomy, in 
particular with SIRTF. 

2. MID INFRARED SURVEYS 

2.1. Specificities 

The mid-infrared surveys have been conducted with ISO- 
CAM (Cesarsky et ah, 1996), mainly at 7 and 15 /xm; a 
few have been performed at 4.5 /xm and 12 /xm (Clements 
et ah, 1999). At 7 /xm using the LW2 filter (e.g. Taniguchi 
et al., 1997; Oliver et ah, 2002), the contamination by 
stars is important and data at other wavelengths (usually 
near infrared) are needed to identify them. The decreasing 
k-correction is less favorable for detecting higher redshift 
galaxies (dash line in Fig. 2, taken from Franceschini et 
al., 2001 ). 

At 15 /xm using the LW3 Alter (e.g Elbaz et al., 1999), 
the stellar contamination is less of a problem and the k- 
correction of galaxies is favorable for detecting sources up 
to redshifts around 1.4 (Fig. 2). Thus, a significant part of 
the extragalactic results are based on 15 /xm observations. 

2.2. Source Counts 

Source counts at 15 /xm show an impressive consistency 
over four decades of flux between surveys of different depths 
(e.g. Elbaz et ah, 1999, Gruppioni et al, 2002). Data come 
from large and shallow surveys like ELAIS (e.g. Serjeant 
et al, 2000), to narrower and deeper surveys (e.g. Ser¬ 
jeant et al, 1997, Aussel et al, 1999; Altieri et al, 1999). 
The number counts (summarized in Fig. 3 by Elbaz et al.. 
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Figure 3. Differential source counts (normalized to Eu¬ 
clidean) at 15 fim from Elbaz et al (1999). The shadowed 
area represents a scenario with no evolution. Symbols rep¬ 
resent different surveys. 

1999) are in excess by a factor of 10, at the 0.5 mJy level, 
compared to non evolution scenarios. The slope is steep 
(a = 3.0±0.1 in the differential counts, for S^, > 0.4 mJy), 
and a turnover appears at S^, = 0.4 mJy. At fainter levels, 
the convergence (decrease) is fast. 

2.3. Nature of the Galaxies 

With a median redshift of 2 ~ 0.8 and most of the sources 
with 0.5 < z < 2 (Aussel et al., 1999), the MIR sources 
have a mean luminosity of 6 x lO^^L© (Elbaz et ah, 2002). 
Most of them are experiencing intense stellar formation 
of about 100 M 0 yr“^, which appears to be uncorrelated 
with the faint blue galaxy population dominating the op¬ 
tical counts at z~0.7 (Ellis 1997; Elbaz et al. 1999). Fig. 4 
(from Elbaz et ah, 2002) represents their HDF-N sample, 
and summarizes the properties of MIR galaxies: redshift 
distribution, luminosity, and SFR (for 80% of the non- 
AGN galaxies). 

2.4. Galaxy Evolution 

Elbaz et al (2002) present a summary of the galaxy evo¬ 
lution at 15 /xm and beyond. First, they show that 60% 
of the GIB at 15 /xm is created by LIRGs {L > IO^^Lq). 
Second, they show that the comoving luminosity density 
at 15 /xm was about 55 times larger at z ~ 1 than today. 
Third, still using some assumption about the SED of the 
galaxies, they conclude that the comoving density of IR 
luminosity radiated by dusty starbursts (which is directly 
related to the star formation) was about 70 ± 35 times 
larger at z ~ 1 than today. At intermediate redshifts, Flo¬ 
res et al (1999) also determined the star formation rate. 
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Figure J;.. Infrared luminosity (left axis) and star formation 
rate (right axis, for starbursts only) of 15 fum S > 0.1 
mJy galaxies in the HDF-N from Elbaz et al (2002). The 
dash curve represents the sensitivity limit. This figure also 
illustrates the redshift distribution of the sources, with a 
median redshift of z = 0.8 (Aussel et al., 1999). 

Finally, they show that about half of the CIB at 140 /zm 
is produced by LIRGs and about one third by ULIRGs. 
Thus MIR galaxies observed by IS OGAM have resolved 
about 75% of the CIB at 140 pLm. 

3. FAR INFRARED SURVEYS 

3.1. Specificities 

The far-infrared surveys have been conducted with ISO¬ 
PHOT (Lemke et ah, 1996), mainly at 90 and 170 a 
few have been performed at 60, 120, 150 and 180 pLm (Ju- 
vela et ah, 2000; Linden-Voernle et ah, 2000). In the 50 to 
100 p,m range, the ClOO Ge:Ga camera is sensitive to the 
peak of rest-frame emission from obscured star formation. 
The decreasing k-correction, in addition to a challenging 
data processing, led to probe mainly the local universe. 
With 46 arcseconds per pixel and a FWHM of about the 
same size, the angular resolution at 100 pLm is significant 
improved since IRAS. 

In the 100 to 200 pirn range, the C200 stressed Ge:Ga 
camera is sensitive to the very cool local galaxies as well 
as higher redshift starburst galaxies, thanks to an advan¬ 
tageous k-correction (e.g. Guiderdoni et ah, 1997) and the 
behavior of the camera (e.g. Lagache & Dole, 2001). With 
92 arcseconds per pixel and a FWHM of about the same 
size, the angular resolution at 170 pLm is an issue for source 
identifications. 

3.2. Source Counts 

Source counts at 170 /rm (e.g. Kawara et ah, 1998; Puget 
et ah, 1999; Dole et ah, 2001) exhibit a steep slope of 


Figure 5. Distribution of IR luminosities from Patris et al 
(2002). Solid line: FIRBACK FSM 170 pm sample (S > 
200 mJy); Dash: IRAS 100 pm sample from Berlin et al 
(1997) with no < feo < 200 mJy. 

a = 3.3 ± 0.6 between 180 and 500 mJy (lower left panel 
on Fig. 7) and, like in the MIR range, show sources in 
excess by a factor of 10 compared with no evolution sce¬ 
nario. This strong evolution was unexpected and quite dif¬ 
ficult to reproduce with models. Matsuhara et al (2002) 
extended the source count analysis to fainter fluxes and 
still detected a strong evolution. 

At 60 and 90 pm, the situation is less clear, but pre¬ 
liminary source counts (Kawara et ah, 1998; Efstathiou 
et ah, 2000; Linden-Voernle et ah, 2000) are compatible 
with no evolution scenarios on their bright end, and be¬ 
gin to show evolutionary effects on their bright end. New 
processing techniques (see Sect. 5) will certainly allow to 
go deeper. 

3.3. Nature of the Galaxies 

Determining the nature of the FIR galaxies has been a 
longer process than in the MIR, mainly because of the 
difficulty to find the shorter wavelength counterparts in a 
large beam. Various techniques have been used to over¬ 
come this problem, one of the most successful being the 
identification using 20cm radio data (e.g. Ciliegi et al., 
1999). Another technique is the FIR multiwavelength ap¬ 
proach (Juvela et al., 2000) that helps constraining the 
position and the SED; it also helps to separate the cirrus 
structures from the extragalactic sources. A variation is to 
use ISOCAM and ISOPHOT data, like the ELAIS Survey 
(Oliver et al., 2000; Serjeant et al., 2000,2001). Finally, the 
Serendipity Survey (Stickel et al., 1998, 2000), by cover¬ 
ing large and shallow areas, allows to detect many bright 
objects easier to follow-up or already known. 

FIR ISO galaxies can be sorted schematically into two 
populations. First, the low redshift sources, typically z < 
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Figure 6. Fluctuations of the CIB in a power speetrum 
analysis of the FIRBACK/ELAIS N2 field at 170 microns 
by Puget & Lagaehe (2000). Observed power speetrum: di¬ 
amond; straight eontinuous line: the best fit cirrus power 
spectrum; dash line: cirrus power spectrum deduced from 
Miville-Deschenes et al (2002); continuous curve: detector 
noise. 


0.3 (e.g Serjeant et al., 2001; Patris et al., 2002, Kakazu et 
al., 2002), have moderate IR luminosities, below KA^Lq, 
and are cold (Stickel et al., 2000). Second, sources at 
higher redshift, z ^ 0.3 (Patris et al., 2002, see Fig. 5) 
and beyond, z ~ 0.9 (Chapman et al., 2002) are more 
luminous, typically L > lO^^L©, and appear to be cold. 
Serjeant et al. (2001) derived the Luminosity Function at 
90 pLm, and started to detect an evolution compared to 
the local IRAS 100 pm sample. 

3.4. Fluctuations of the CIB 

Sources below the detection limit of a survey create fluctu¬ 
ations. If the detection limit does not allow to resolve the 
sources dominating the CIB intensity, which is the case in 
the FIR with ISO, characterizing these fluctuations gives 
very interesting information on the spatial correlations of 
these unresolved sources of cosmological significance. An 
example of the modeled redshift distribution of the un¬ 
resolved sources at 170 pm can be found in Fig. 12 of 
Lagaehe et al (2002); the sources dominating the CIB fluc¬ 
tuations have a redshift distribution peaking at z ~ 0.9. 
After the pioneering work of Herbstmeierer al (1998) with 
ISOPHOT, Lagaehe & Puget (2002) discovered them at 
170 pm in the FIRBACK data, followed by other works at 
170 and 90 pm (Matsuhara et al., 2000; Puget & Lagaehe, 
2000; Kiss et al, 2001). Fig. 6 shows the CIB fluctuations 
in the FN2 field by Puget & Lagaehe, (2000), at wavenum¬ 
bers 0.07 < k < 0.4 aremin”^. 


4. MODELS 

ISO provided high quality data, taken in the frame of the 
cosmological survey programs, as well as nearby galax¬ 
ies or serendipitous programs. The first observable to be 
widely used as a constraint on the models is the source 
counts at 15 and 170 pm, but also at 7 and 90 pm. The 
redshift distributions, mainly at 15 pm (and soon at 170 
pm) put additional constraints, as well as a better knowl¬ 
edge of the nature of the galaxies: AGN vs Starburst, lu¬ 
minosities, star formation rate. The global star formation 
rate, as derived in part from ISOCAM data, adds other 
constraints to the models. Finally, a few models use ex¬ 
plicitly the constraint of the level of the CIB fluctuations 
in the far infrared. 

These observables and constraints induced a burst in 
the development of new models in the late nineties, among 
which: Roche & Bales (1999), Tan et al. (1999), Devriendt 
& Guiderdoni (2000), Dole et al. (2000), Wang (2000), 
Chary & Elbaz (2001), Franceschini et al. (2001), Malkan 
& Stecker (2001), Pearson (2001), Rowan-Robinson (2001), 
Takeuchi (2001), Xu et al. (2001), Balland et al. (2002), 
Lagaehe et al. (2002), Totani & Takeuchi (2002), Wang 
(2002). This development was also made possible because 
of the better knowledge of the galaxies SED in the infrared 
(e.g. Dale et ah, 2000; Helou et ah, 2000; Lutz et al., 2000; 
Stickel et ah, 2000, Tuffs & Popescu, 2002), with the avail¬ 
ability of the SED of the CIB (Gispert et ah, 2000; Hauser 
& Dwek, 2001), and of complementary data sets (X-rays, 
optical, NIR, submm, radio). 

Even if there is no unique scenario reproducing all the 
observables, the models help now to have a more coherent 
picture of galaxy evolution. It is beyond the scope of this 
paper to review the models and their predictions, we thus 
give here only two examples. Fig. 7 shows the observed 
source counts at 15, 60, 170 and 850 pm and the fits by 
Lagaehe et al (2002). Using other kind but similar evo¬ 
lutions (pure density and pure luminosity evolutions in 
addition to the density-1-luminosity evolution), Chary & 
Elbaz (2001) predict the star formation rate and compare 
it with data and Xu et al (2001) (Fig. 8). Combined with 
the different approach of Gispert et al (2000), the SFR is 
well constrained in the redshift range 1 < z < 3. 

5. POTENTIAL OF ISO DATA 

Even if the ISO data led to an impressive scientific re¬ 
turn with a high efficiency (about 62% of the data are 
published) in the field of the cosmological evolution of the 
galaxies, many other results are still to come, based on 
the amount of unpublished data available in the archive. 

In addition, new techniques of data reduction are now 
available, thanks to a better understanding of detector 
behaviors, both for the Silicon Array (CAM)(Coulais & 
Abergel, 2000; Miville-Deschenes et ah, 2000; Lari et al., 
2001, 2002; Vaccari, 2002) and the Germanium detectors 
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Figure!. Observed source counts at 15, 60, 170 and 850 
fim with the model of Lag ache et al (2002). Symbols repre¬ 
sent data (CAM at 15 pm, IRAS at 60 pm, PHOT at 170 
pm and SCUBA at 850 pm. Solid line: model prediction. 
Dash: LIRC/starburst component. Dot-dash: normal/cold 
galaxy component. 


(PHOT) (Coulais et al., 2000). Processing unpublished 
data and/or reprocessing published data allows us to reach 
better sensitivity levels with a higher confidence (e.g. for 
ISOCAM HDF-N: Serjeant et al., 1997; Aussel et al., 1999; 
Desert et al., 1999. E.g. for ISOPHOT Lockman Hole: 
Kawara at al., 1998; Rodighiero et al., 2002). 

Concerning the CIB and the cirrus foreground, a bril¬ 
liant example of what can be done from the archive is the 
work of Kiss et al (2001,2002), detailed in this conference. 
Furthermore, the ISO data may induce discoveries in other 
data sets. A recent example concerns the fluctuations of 
the CIB in the far infrared. Lagache & Puget (2000) and 
Puget & Lagache (2000) measured the level of the fluc¬ 
tuations of the CIB at YJDpm in the FIRBACK fields. 
Given the SED of the CIB, it was possible to predict the 
fluctuation level at 60 and 100 pm (e.g. Lagache, Dole, 
Puget, 2002), and show that this might be detectable in 
the IRAS data. As a consequence, Miville-Deschenes et 
al. (2002) recently reported the first detection of the CIB 
fluctuations at 60 and 100 pm in the IRAS data. 


6. THE NEXT STEPS: SIRTF AND BEYOND 


SIRTF, the Space Infrared Telescope Facility, to be laun¬ 
ched in early 2003, is expected to improve our view of the 
IR universe over the next five years, thanks to the recent 
technological developments in IR detector arrays. A num¬ 
ber of multiwavelength cosmological surveys are scheduled 
with MIPS (Multiband Imaging Photometer for SIRTF) 
at 24, 70 and 160 pm (Dole et al., 2003[]) and IRAC (In¬ 
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Figures. Star Formation Rate from Chary & Elbaz (2001). 
Upper panel: min and max range from their model, and 
observed UV/opt data; dots represent Xu et al (2001) 
model. Lower panel: 3 different evolution scenarios from 
their model and data corrected for extinction. Line: pure 
luminosity; upper dash: pure density; lower dash: luminos¬ 
ity-!-density. 


frared Array Camera) at 3.6, 4.5, 5.8 and 8.0 pm. In the 
frame of the Guaranteed Time (GTO) or Legacy (GOODS 
and SWIRE) Surveys, about 80 Sq. Deg. will be observed 
at different depths. Probing the properties of the galaxies 
with samples of statistical significance up to redshifts of 
2.5 or more will be possible. Fig. 9 (from Dole et al, 2003) 
shows predictions of the redshift distribution at 24 pm for 
MIPS surveys. 

The next steps after ISO and SIRTF for understanding 
the infrared and submillimeter Universe from space will 
be ISAS’s ASTRO-F in 2004, ESA’s Herschel & Planck 
(around 2007), and NASA’s NGST (around 2010-I-), in 
addition to the ground-based ALMA (from 2005). These 
telescopes are the promise for a lot of exciting science to 
be done (at least) in the next 15 years ! 
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Figure 9. Predicted redshift distribution of sources observed 
at 24 tJ-m with MIPS on SIRTF, from Dole et al (2003), 
for various 5 (t surveys depths. Line: GTO shallow sur¬ 
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